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Abstract 

For bouncing cosmologies such as the ekpyrotic/cyclic scenarios we show that it is possible to make 
predictions for density perturbations which are independent of the details of the bouncing phase. 
This can be achieved, as in inflationary cosmology, thanks to the existence of a dynamical attractor, 
which makes local observables equal to the unperturbed solution up to exponentially small terms. 
Assuming that the physics of the bounce is not extremely sensitive to these corrections, perturbations 
can be evolved even at non-linear level. The resulting spectrum is not scale invariant and thus 
incompatible with experimental data. This can be explicitly shown in synchronous gauge where, 
contrary to what happens in the commonly used Newtonian gauge, all perturbations remain small 
going towards the bounce and the existence of the attractor is manifest. 



1 Introduction 

From cosmological observations we know that the current Universe is to a good approximation flat, 
homogeneous and isotropic on large scales. It is well known that in the standard Big Bang cosmology 
this requires an enormous amount of fine tuning on the initial state emerging from the Planck 
era. This problem is usually solved postulating a period of inflation which acts like a dynamical 
attractor for the cosmological evolution. From a generic initial condition sufficiently close to a flat, 
homogeneous and isotropic state the Universe evolves towards this symmetric configuration. Another 
logical possibility that has been explored is that a dynamical attractor is present in a contracting 
phase which precedes the present expansion. After this phase, the Universe would reach a state of 
high curvature and then bounce and start expanding. Although string- inspired, the ekpyrotic/cyclic 
scenarios 12] are fully described in the contracting phase well before the bounce by a 4d effective 
field theory. The only light degree of freedom (besides the graviton) is a scalar field which moves along 
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a steep negative potential and drives the cosmic contraction . Even in the presence of perturbations, 
if the potential is steep enough, inhomogeneities and anisotropies become less and less relevant going 
towards the bounce and the cosmic evolution gets closer and closer to the unperturbed one. 

The crucial question is whether this class of models can give rise to an approximately scale 
invariant spectrum of adiabatic density perturbations as clearly required by experiments. The answer 
to this question is not straightforward because there is no explicit model of a bouncing phase; it is 
just assumed that at sufficiently high energy UV corrections will stop the contraction and lead to 
an expanding phase. It is therefore important to understand if one can make predictions which are 
robust, i.e. independent of the details of the unknown bouncing phase. In inflationary cosmology 
we face a closely related issue because perturbations, that are created during inflation, are observed 
much later and there are cosmological phases in between, like reheating, whose details are completely 
unknown. It is well understood that in models with a single fluctuating field all predictions are 
independent of the details of the unknown cosmological phases. There are two reasons why this is 
true. First of all we are interested in modes which are well outside the horizon during the unknown 
regimes of the cosmological evolution. This would not be enough however. The situation is further 
simplified by the absence of isocurvature perturbations: every observer will eventually go through the 
same cosmological history. Following the dynamical attractor, after a mode goes outside the horizon, 
the Universe locally approaches the unperturbed solution up to exponentially small terms. This 
"parallel Universes" approach ^4 allows to follow perturbations through the unknown cosmological 
regimes at any order in perturbation theory. On the contrary if we add another field into the game, 
separate regions of the Universe will be characterized by different values of this new field and will thus 
undergo a different evolution, so that the final state depends on the details of the whole cosmological 
history. The fate of the isocurvature component can be very different: it can be washed away by 
thermal equilibrium or, on the contrary, become the leading source of adiabatic perturbations as in 
the curvaton and variable decay width scenarios. 

The purpose of this paper is to understand whether the same kind of arguments used for single 
field inflation can be applied to evolve density perturbations through the unknown bouncing phase 
of the ekpyrotic/cyclic scenario. 

In Newtonian gauge the most generic scalar perturbation is (at linear order) a superposition of the 
two independent solutions for the Newtonian potential It has been pointed out in the context of 
the ekpyrotic/cyclic scenario [J] that, in the limit of a very steep potential, one of these two solutions 
has an approximately scale invariant spectrum in the contracting phase. Unfortunately the most 
naive evolution from the contracting to the expanding phase suggests that this scale invariant mode 
is orthogonal to the "growing mode" in the expanding phase, the one relevant for observations. 
However it has been claimed that a generic mixing of the two solutions at the bounce would induce 
an observationally viable spectrum. In this paper (section OJ) we show that this generality can not 
be defended in any way. For example if we change gauge and describe the most generic perturbation 
in terms of the variable £ (proportional to the curvature of comoving surfaces), none of the two 
independent modes is scale invariant. This tells us that for a contracting cosmology the concept 
of scale invariance of a particular variable has no physical meaning before one specifies how this 

the pre-Big Bang scenario [3] the (Einstein frame) contraction is driven by the kinetic energy of the 
dilaton, which has no potential. As we will discuss this can be seen as a limiting case of the ekpyrotic/cyclic 
models and the conclusions wc will draw arc different. 
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variable is related to what we will eventually observe. 

As we discussed, in the presence of a dynamical attractor in the contracting phase, close to 
the bounce the cosmological solution is locally (i.e. on regions smaller than the Hubble radius) 
homogeneous, isotropic and flat, and the scalar follows its unperturbed evolution, up to exponentially 
small terms. This leads us to make a simple assumption (section HJ), which is implicitly made also 
in the case of an inflationary Universe when evolving perturbations through unknown phases. If 
the physics of the bounce is not tremendously sensitive to these exponentially small corrections, 
every observer will go through the bounce in the same way, following the unperturbed evolution. 
This allows to predict, not only at linear level but at any order in the perturbation, the statistical 
properties of the fluctuations in the expanding phase. Unfortunately the conclusion is that the 
spectrum of density fluctuations is not scale invariant and thus it is ruled out by experiments. This 
can be explicitly shown in synchronous gauge where, contrary to what happens in the commonly 
used Newtonian gauge, all perturbations remain small going towards the bounce and the existence 
of the attractor is manifest. If our assumptions are not satisfied, which is a logical possibility, then 
all predictions, like the scale invariance of the 2-point function or the level of non-gaussianity of the 
perturbations, strongly depend on the details of the bounce and no robust prediction can be made. 

We stress that we are not giving any prescription to match perturbations across the bounce 
(for example we are not claiming that ( must be continuous [HI IH1 E] j although the final result 
for density perturbations in the expanding phase is the same). Given a generic initial condition 
before the bounce there is no way to evolve it through the unknown phase without specifying all 
the details. The point here is that we do not have a generic initial condition but a very special 
one, because of the presence of the dynamical attractor. Locally the solution is exponentially close 
to the unperturbed solution and this is a powerful simplification. Again the situation is similar 
to inflationary cosmology; nobody knows how to evolve a generic initial state through reheating 
because we do not know how reheating looks like, but there is no need for this as inflation leads to a 
very specific state, locally (but not on large scales) indistinguishable from the unperturbed solution. 



In this section we want to approach the problem by a (critical) review of the basic heuristic argument 
for a scale-invariant spectrum of density perturbations in a bouncing cosmology [51 1101 ITT] . This 
will show that there are qualitative differences with respect to an inflationary scenario that must be 
taken into account. For the time being we neglect how perturbations will evolve through the bounce. 
The argument is an estimate of the fluctuations of the scalar field 4>, which is driving the contraction 
of the Universe towards the bounce. As the field is pulled by a steep negative potential, it seems 
natural to assume that gravity is negligible in the dynamics of the fluctuations of (j). Extending the 
discussion of refs. [El El, we will argue that this is not the case. 

To simplify the algebra it is useful to approximate the potential over some range of (j> with the 
exponential form [2J 



In the ekpyrotic/cyclic scenario the potential is very steep, i.e. p <C 1. The Friedmann equations 



2 Intuitive argument and its limitations 
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and the equation of motion for an homogeneous (ft configuration, 



4> + 3H<p + V'{(l)) = , (2) 

are exactly solved by the background 

a(t) = (t/t f, H(t) = Z, Mt) = M PV % log I , (3) 

where to = — yPO- ~ ' Mp/t/Vq. The scalar <f> is moving from +oo to — oo and t is negative and 
running towards which corresponds to the unknown bouncing phase. This solution has a constant, 
large pressure-to-energy density ratio, w = (2/3p) — 1^1. 

It can be shown that for p < 1/3 this solution is a dynamical attractor, in the sense that the 
system rapidly approaches it starting from a generic homogeneous and isotropic initial condition 
[T2*] . Roughly speaking, in the presence of a very steep potential the kinetic energy of the scalar 
becomes bigger and bigger in the evolution, so that the asymptotic state does not depend on the 
initial kinetic energy. In the following we will explicitly see that the same conclusion remains true 
in the presence of (sufficiently small) departures from isotropy and homogeneity. As discussed, this 
property will turn out to be crucial for our argument. 

Since in the limit p <C 1 the scale factor changes very slowly with time, one could be tempted to 
say that the gravitational backreaction is small and study perturbations in the scalar cj> on the fixed 
gravity background. The equation of motion on a fixed FRW metric would be 

<% + 3#<% + + V"{4>)\ <% = , (4) 

where k is the comoving wavevector. For the background above we get 

3p- (k 2 2(1 -3p) 



<% + f % +{-2- ) <% = • (5) 

It is evident that the contraction of the Universe gives a subleading contribution for p <C 1 • Moreover 
for small t the gradient term can be neglected so that the mode freezes. The equation reduces to 
— j?d<fih = 0) which gives the two behaviors 5(j)? ~ i -1 and t 2 . Therefore the leading solution for 
small t, properly normalized to match the usual Minkowski limit for large negative t, is 

This implies a scale-invariant spectrum in position space as (5(J)(x) 2 ) ~ J dlogk k^\5(j)^\ 2 . 

It is however easy to realize that the argument above is not physically motivated, even at the 
heuristic level. In fact there is no limit in which the fluctuations of <j) can be decoupled from the 
metric perturbations since <f> is what drives the evolution of the Universe and its potential is steep. 
Taking the naive decoupling limit Mp — > oo one obviously ends up with a trivial background. The 
only correct decoupling limit would be to have a flat potential for (ft; in this case eq. (jlj is correct 
because fluctuations of (j) do not gravitate except for gradient terms. This is the reason why in 
slow-roll inflation the intuitive argument above gives a correct estimate of the density perturbations 
up to slow-roll corrections. 
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In our case there is no sense in which the dynamics of perturbations is dominated by the potential 
term: gravity is a crucial ingredient and must be taken into account. Let us clarify this statement in 
a gauge which is quite close to the intuition of keeping 5(f) as the dynamical variable characterizing 
scalar perturbations, 



Unless otherwise specified we concentrate to scalar modes, setting tensor modes to zero. In this 
gauge gravitational effects are suppressed in inflation by slow-roll parameters, so that it explicitly 
realizes the decoupling limit discussed above. 

The other components of the metric, goo and goi, are Lagrange multipliers and as such are 
determined as a function of 5(f) through the constraint Einstein equations. Substituting them into 
the action for gravity and the scalar field and specializing to our background solution we get a 
Lagrangian for 5(f) of the form |13j (see Appendix lAl for the derivation of this result) 



This is the Lagrangian for a scalar in the given background but without any potential term, so that 
from this point of view it is difficult to argue that the potential is driving the dynamics of the 
perturbations. Of course the equation of motion deriving from this Lagrangian is exactly eq. (|2j) 
without the mass term. The effect of metric perturbations is so important that exactly cancels 
out the potential. Note that the very simple Lagrangian above is exact only in a scaling solution 
a(t) tx t p , like the collapsing background described above (p <C 1) or inflation with an exponential 
potential (p> 1). In both cases the action will receive small corrections proportional to the variation 
of the equation of state in one Hubble time, i.e. subleading respectively in the fast-roll or slow-roll 
expansion. 

From the action above, the properly normalized solutions for small t are now 5(f)^ ~ and 
5<p^ ~ \fk t. None of these implies a scale-invariant spectrum. 

To get further insight into the problem, in the next section we will perform a detailed analysis 
of scalar perturbations in two different gauges. 



3 Standard gauges and the apparent rapid growth of 
perturbations 



The study of scalar perturbations can be performed in different gauges. In each one, after satisfying 
the constraint equations, the perturbation is parametrized by a single scalar function, which satisfies, 
at linear order, a second order linear differential equation. The most generic scalar fluctuation will 
therefore be described by a linear combination of the two independent solutions of this equation. 
Usually one of the solutions dominates over the other at late times ( "growing" and "decaying" mode). 
However in some cases the same physical perturbation is described in one gauge as growing mode 
and in another as decaying (we will see an example of this in the following) . Neglecting the decaying 
mode is therefore quite misleading, especially in a contracting background like the one at hand, and 
thus will be avoided. 

To make contact with the literature, in this section we are going to perform the calculation 
in two commonly used gauges. Doing so we will stress some subtle point which has been usually 





(7) 




(8) 
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overlooked. However it will turn out that in both gauges perturbation theory breaks down for t — > 
in the background under study. In the next section we will therefore move to synchronous gauge in 
which the perturbativity of the calculation remains manifest. 

The analysis of perturbations in inflation is greatly simplified by using the variable £ introduced 
in ref. because it is conserved outside the horizon under simple assumptions on the stress energy 
tensor. We therefore start our analysis in the ("-gauge defined by 

cft(x, t) = Mt) , 9ij = (1 + t)) a 2 (t)<% . (9) 

In this gauge the scalar field is unperturbed so that the relation with the gauge introduced in the 
previous section is simply a time reparameterization 

f-t-^ =* c = -5V (io) 

In our scaling background solution 5 (ft and ( are simply related by a constant factor 

e = _,/£JL$0 (ii) 

s V 2 M P Y y ' 

so that the Lagrangian for Q is again very simple 

M 2 f 

S = £ d 3 xdtV=g(d^() 2 . (12) 

p J 

The other components of the metric are fixed as functions of £ through the constraint equations. At 
linear order they are (see Appendix ^] for details) 

<7nn - I = — 1 — — C (13) 

where, as usual, 1/d 2 is defined in Fourier transform. Note that £ is related to the Ricci scalar of 
the spatial metric induced on comoving surfaces (i.e. spatial hypersurfaces of constant (ft) simply by 

(3) R = 4& > (15) 

as in this gauge comoving surfaces are surfaces of constant t. 

The equation of motion for £ deriving from its action above is 

Q n Z-.2 

4 + T ( E + W# ( "' = 1 ' (16> 

which can be integrated in terms of J and Y Bessel functions with argument kr, where r is the 
conformal time defined as usual by dt = a(t)dr. The full solution is 

^(r) = -V^t? (Aj-Ji^i-kr) + B^Y^i-kr)) , (17) 

CL \ 2—2p 2~2p / 

where and are integration constants, to be determined below. It has been claimed that the 
variable £ is not suitable for discussing density perturbations in a contracting phase, because it shows 
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no "classical instability" and therefore is not "amplified" [HUH]. The meaning of this is unclear since 
the two modes for £ above parametrize the most generic scalar perturbation of the system. 

Following the standard procedure TE\ , we then quantize the system by writing the quantum field 
operator as CjfCr) = Cfc( r ) a g + C?( T )^ r> where C^( r ) is the above classical solution, and the creation- 
annihilation operators obey the standard commutation rule [a^, at] = (27r) 3 5 3 (fc — k'). Since for 
large negative r the Universe reduces to Minkowski spacetime, we impose that the classical solution 
Ct( T ) asymptotes to the standard Minkowski wavefunction oc -hr^~ lkT of a massless scalar field. 



Also, we assume that the initial state is the usual Minkowski vacuum. The proper normalization 
can be read from the action eq. ()12j) written in conformal time, 

Cr(r) — ■ -^=e~ ikr for r -oo . (18) 

k ' V2M P a ^2k V ' 

This fixes the linear combination of the two independent solutions to be 

Cg(r) = - 1-^ ■^V^{j,_ 2R {-kr)+iY 1 ^ E (-kr)) . (19) 
At small t we recover the two behaviors that solve eq. (|16j) neglecting gradient terms 

C * ~ W P {lk ~ l+ ' % + ^"^ 1_3P ^) , (20) 

where we expanded the exponents up to linear order in p. The first term comes from the Y function 
while the second comes, for generic p, both from Y and J. Even though £ remains finite for t — > 
the off-diagonal term goi in eq. (|14j) diverges like t~ p , so that linear theory cannot be justified in this 
gauge for small t. 

The late time spectrum of C is 

(4(i)4-,(t)> = {2vf5\k + k') KfeWI 2 , (21) 

with 

|C^)| 2 ~ ^ik-^tf + k^-^tf) , (22) 

which does not contain any scale-invariant component. This remains true even if one includes all 
further terms in the kr expansion of the complete solution eq. (|19|) as they only add positive powers 
of k. This result can be written in a more transparent form as a function of the physical wavevector 
^-ph = kj a, 

<C(*) 2 ) ~ f dlog k ph pf^ [(k ph /H) 2+2p + (k ph /H) 4 ~ 2p ] . (23) 

We stress again that the complete solution for £ describes the most generic scalar fluctuation 
of the system. But, as discussed in the Introduction, a scale invariant spectrum is obtained for the 
Newtonian potential <& so that the situation is at first puzzling. 

Before analyzing the relation between the results obtained for ( and <£, it is important to appre- 
ciate that the second term in eq. 1(2*0]) cannot be simply disregarded as "decaying" . Although it drops 
off in time like i 1_3p , thus giving a negligible contribution to the intrinsic curvature of comoving 
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surfaces (see eq. (|15p). it dominates the perturbation in the extrinsic curvature of comoving surfaces. 
In fact at linear order the extrinsic curvature is 



K*j = y\g kj - d k90j - d j9h0 ) = (H + 06} - -^C - -^C • (24) 

J p o z p a z 

The leading contribution comes from the term £c ~ i 1_3p in eq. (|20[). which gives a perturbation 
in K l j diverging as t~ 3p towards the bounce. Note that this is however a small perturbation with 
respect to the unperturbed extrinsic curvature which diverges faster, as H ~ t . We see that the 
two terms in eq. (1201) are both physically relevant as they dominate respectively the intrinsic and 
extrinsic curvature. Further terms in the small t expansion of eq. (|19|) are instead subleading in 
their contribution to physical curvatures. 

Let us now see how these perturbations are reinterpreted in the commonly employed Newtonian 
gauge. This gauge is defined by 

ds 2 = -(1 + 2<$>)dt 2 + (1 - 2^)a 2 {t)dx 2 . (25) 

The most generic scalar perturbation can be described by the Newtonian potential <F All other 
variables are related to <& through the constraint Einstein equations. In particular \I/ = $ at linear 
order in absence of anisotropic stress, as in our case. The equation of motion for $ in our background 
is (see for example ref. [3]) 

*» +? 7 £ *« + ci7S55»i- - (26> 

Following the usual procedure as done for Q we get the properly normalized solution 

^ k {r) = l -^^.-^^=(j 1± ^(-kT)+ l Y 1±]L {-kT)) , (27) 
a V2/c 2Mp 1 - p yj-kr V 2-2 P 2-2 P / 

with late time behavior 

~ A-{k- z/2 - p t- 1 ~ p + k-^ 2+p t p ) . (28) 

We see that the first term, which diverges approaching the bounce (thus casting doubt on the linear 
approximation), has an approximately scale invariant dependence on k as first pointed out in ref. 
( 2 ). But how is it possible that here we have an almost scale invariant spectrum, while there is no 
sign of scale-invariance in £? After all, we said that the most generic perturbation can be described 
using either £ or <F It is straightforward to verify that the scale invariant term in the expression 
for dominates the extrinsic curvature of comoving surfaces. In fact, from the relation between 8<p 
and <3? in this gauge [5] 

2M 2 

6(f) = -^(<$> + H$>) , (29) 

we obtain the contribution of the first term in eq. (|28|) to the extrinsic curvature of the 5<j) = 
surfaces 

Kij-H-^Qzoct-V 2 , (30) 

2 Note that even though this seems to validate the naive argument discussed in the previous section, the 
time dependence of the second solution, $ = const, is not a solution of eq. (J^J. The same happens for 8<j> in 
Newtonian gauge which has the same time dependence as $ (see eq. J23))- 
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which is the same t and k dependence we obtained in the previous gauge, eq. (EU) ( 3 ). Notice 
that the different dependence on k of the two variables £ and $ corresponds to the fact that their 
relation with physical quantities like curvatures can involve a different number of derivatives. This 
tells us that the scale invariance of a scalar fluctuation has no well defined meaning before specifying 
which variable is relevant for observation. The scale-invariant term in eq. (|28j) gives a vanishing 
contribution to the intrinsic curvature eq. (|15|). This is clear from the explicit relation between £ 
and $ 



that is why there is no sign of scale- invariance in (. The fact that £ vanishes for the first term in 
eq. (|2*5)l could suggest that we are losing it in going to the variable £, but as we saw its geometric 
effect in perturbing the extrinsic curvature is completely captured by £. This is also a nice example of 
how misleading the description of a solution as decaying or growing can be; the behavior <I> ~ i~ 1-p 
would be considered as growing in Newtonian gauge, but we saw that the same physical fluctuation 
is "decaying" in ( gauge, ( ~ t 1_3p . 



4 Synchronous gauge: making the dynamical attractor 
manifest 



We are now going to discuss the relation between the scalar perturbations produced in the contracting 
phase and what is observed today. The behavior at late times of scalar fluctuations in an expanding 
phase can be obtained from our eqs. (|2U|) and (|28[) (with generic coefficients), remembering that 
now t is positive and going towards +00. Note that now the time independent piece in both <3? and 
£ dominates each variable, with the time dependent piece going to zero (now p = 2/3 in matter 
dominance and p = 1/2 in radiation dominance). 

In the present expanding phase all the ambiguity described above about the different intuition in 
the two gauges and the distinction between growing and decaying modes is absent. This is because 
once the fluctuation reenters the horizon all physical quantities depend (apart from exponentially 
small corrections) only on the £ ~ <£> ~ const mode; the concept of growing and decaying modes has 
thus a clear physical meaning. To have a viable model the constant mode must have an approximately 
scale invariant spectrum. Even though the time independent term of in eq. 1)28(1 does not have a 
scale invariant spectrum in the contracting phase, it has been argued that the bounce can induce a 
generic "mixing" between the two solutions, so that the scale invariant spectrum of the $ ~ t~ l ~' p 
piece before the bounce is inherited also by the $ = const solution in the expanding phase. It is 
however clear from the discussion above that it is difficult to defend this point of view, as in the 
C-gauge there is no term with a scale invariant spectrum. It is therefore crucial to identify the 
relevant variable in which this mixing, if any, happens. 

3 To be precise only the anisotropic piece of the extrinsic curvature in eq. (|24|l is recovered in this way. 
To get the isotropic component proportional to £ we have to keep further terms (~ i 1_3p ) in the small t 
expansion of $ in addition to those of eq. I|28[> . This again tells us that it is not obvious to get the physical 
importance of a term by simply looking at its contribution to <I>; a term diverging as i _1_p and one decaying 
as t 1_3p give comparable extrinsic curvatures! 




(31) 
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The first issue one should address in studying the evolution of the perturbations towards the 
bounce is the validity of perturbation theory. In Newtonian gauge both $ and 5(f) blow up as 
in the £-gauge, even though £ remains finite all the way to the bounce, the goi components of 
the metric eq. (|14|) diverge as t~ p . These facts indicate that a fully non-linear description might be 
required to approach the bounce in these gauges. Note however that as t — ► the physical curvatures 
induced by the perturbations are smaller and smaller with respect to H, which sets the scale of the 
unperturbed curvatures. Similarly to the inflationary scenario, the Universe becomes closer and 
closer to the unperturbed solution, as implied by the existence of an attractor. Thus, it is obviously 
useful to look for a gauge in which the smallness of perturbations remains manifest approaching the 
bounce. 

This is achieved in synchronous gauge defined by goo = —1, go% = 0. Performing a gauge 
transformation from the £-gauge and specializing to our scaling background we get for the metric 
and the scalar field 4 



ds 1 = -dt 2 + a 2 (t) 



1 + 2C- 2 H J\' - HM C + (C ■ k 2 t 2 ~ 2p ) 



dxUx j (32) 



^ M "f dt 'i. (33) 



t Jo H 

These expressions are explicitly derived in Appendix [Bj In the anisotropic piece of the metric only 
the time dependent part of £ is relevant as a constant term can be set to zero by a spatial gauge 
transformation. The terms in brackets are either constant or decaying in time, making manifest 
that perturbations stay small and we are approaching the unperturbed solution; perturbations are 
becoming locally unobservable. After a mode goes outside the horizon a local observer will still feel 
a residual curvature, with a consequent deviation from the unperturbed FRW solution. A curvature 
term blueshifts as or 2 ~ t~ 2p , so that it becomes rapidly irrelevant with respect to H 2 ~ t 2 ; 
this effect is described by terms going as t 2 ~ 2p in the solution above. In the metric above there 
are also terms decaying as i 1-3p , i.e. slower than curvature terms. They describe the effect of 
anisotropies, which blueshift as a -6 . Their time dependence is compatible with the fact that they 
become irrelevant as long as the unperturbed solution has w > 1 (so that its energy density blueshifts 
faster than a -6 ) as we are assuming. The physical significance of the terms that we have sketched 
will be much more transparent when we will describe the full non-linear solution close to the bounce 
in the next section. Finally also 5(f) goes to zero in this gauge, as i 1_3p ; again this agrees with the 
existence of an attractor for p < 1/3. This shows that the divergence of the Newtonian potential is 
just a gauge artifact; the deviations from the unperturbed solution are getting small. 

We now have a good description of the perturbed metric before the bounce, in which linear 
theory does not break down. We can finally address the issue of how perturbations evolve through 
the bounce to the expanding phase. We are not going to specify anything about the transition 
between the contracting and expanding regimes; our purpose is to understand at what level the 
final predictions are independent of the details of the bounce. The crucial simplification comes from 
the fact that we are interested only in modes with huge wavelength compared to the Hubble radius 



4 The synchronous gauge does not fix completely the reparameterization invariance and one should get rid 
of the remaining gauge modes when writing the equations of motion in this gauge. However the solutions for 
( are here derived in a different gauge. 
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before we enter in the unknown bouncing phase. From the expression above the contribution to any 
local observable O of a given mode with physical wavevector fc pn is suppressed with respect to the 
time-independent part of £ (which we call £ const ) by 

AO C COBSt (^gxV~ 2P (34) 



O * \ H # 

If we assume that the steep potential phase lasts at least 60 e- folds, we obtain that the modes 
of interest have an incredibly small effect for an observer entering the bouncing phase. The ratio 
(fcph/ H) 1 ~ 2p is of order e~ 60 and it describes the fact that the local effect of a mode is very small 
in the presence of the dynamical attractor. Moreover the initial departures from the unperturbed 
solution when the mode leaves the horizon are of order ^fp ■ H cross i ng /Mp (see eq. (J23J)- This gives 
the size of £ const in the equation above. For the modes of interest this gives a further exponential 
suppression e -60 . 

In the concrete proposed models the scaling solution we are describing ceases to be a good 
approximation even before we enter the bouncing phase, because the potential is modified at large 
negative (f). The scaling solution is however valid when all cosmologically relevant modes exit the 
horizon; this ensures that their physical local effect gets exponentially small. 

These considerations lead us to a simple assumption, which makes the predictions of this class of 
models independent of the details of the bounce. We assume that the unknown phase is not sensitive 
to this very tiny effects; every comoving observer goes through exactly the same history (whatever it 
is) apart from exponentially small deviations. With this assumption we can neglect all terms going 
to zero in the metric above 

ds 2 = -dt 2 + a 2 (t)(l + 2( const (x))dx 2 , (35) 

where £ const (x) is the time-independent term of C(x, t). It is important to realize that £ co (x) has an 
exponentially small local effect as, neglecting gradients, it can be locally reabsorbed with a rescaling 
of the spatial coordinates. Now the crucial point is that the form of the metric eq. ()35|) will hold also 
in the expanding phase after the bounce, with the same time-independent rescaling of the spatial 
coordinates £ const (x). Assuming that we have a continuous history through the bounce, any observer 
is going through the same unperturbed history (H(t), 4>(t)) once we have neglected the exponentially 
small corrections. But at each point we have in the metric eq. (|35[) a different normalization of the 
scale factor (l+£ const (x))a(i). This implies that £ const (af) cannot change as it is locally an integration 
constant of the unperturbed solution. The reader could be worried that our conclusions depend on 
the validity of a 4d metric description of the bounce with a continuous solution for the scale factor 
a(t). Actually it is easy to realize that the result holds independently of the unknown physics of the 
bounce, even if it cannot be described in a 4-dimensional field theory language. Given an unperturbed 
solution, every comoving observer will follow it (apart from exponentially small deviations), so that 
the measured red-shift factor between a given comoving time t before the bounce and a given time 
after must be the same for all observers; Q co (x) therefore cannot be modified in the unknown 
phase. 

Our assumption therefore implies that the time independent term of Q is the same before and 
after the bounce. The spectrum of density perturbations at horizon re-entering (Sp/p)* in the present 
expanding phase is thus given by the first term in eq. 1)22(1 . 

<(<Wp)*) ~ (CO oc k- 1+2p . (36) 
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The spectrum is not scale invariant and therefore incompatible with observations. 

We want to stress that our assumption is not trivial. The Hubble constant must change sign 
during the bounce; one can therefore envisage a model in which the bounce is sufficiently "slow" so 
that H remains small for a while and all interesting modes come back into the horizon and start 
oscillating before freezing again From all we said it should be clear that whether or not we 

obtain a scale invariant spectrum does now depend on the unknown physics at the bounce and it is 
not a generic prediction. For example we saw that £ contains no scale invariant term, so that it is 
not possible, without specifying the model, to state that a generic mixing of modes will give a viable 
prediction. 

Another possibility to get around our conclusion without allowing all the modes to come back 
into the horizon is to assume that the physics at the bounce is sensitive to the exponentially small 
corrections we have neglected. The local effect of the perturbations, which has been exponentially 
suppressed during the contraction as a consequence of the attractor behavior, could be amplified 
and "resurrected" . This is physically hard to believe and it would require a certain amount of fine 
tuning to pump up a tremendously small effect and stop at the desired 10~ 5 level before entering the 
non-linear regime. Anyway also this possibility does not generically give a scale invariant spectrum. 

To avoid confusion it is important to underline that our approach does not tell us how to evolve 
a generic scalar perturbation to the expanding phase. A generic fluctuation will be characterized by 
two constants associated with the two independent modes. Different comoving observers will follow 
different histories, for example at the same value of H they will see a different value of the scalar 
field. How the bounce occurs will depend on these differences and to know the final state we must 
obviously specify the details of the bounce. Fortunately the existence of an attractor solution drives 
every observer to see the same history, which is locally indistinguishable from the unperturbed one: 
at the same value of H we have the same value of <j> and <f> apart from exponentially small terms. 
With this crucial simplification the evolution of perturbations through the bounce is fixed by the 
unperturbed solution. 

An approximately scale invariant 2-point function is not the only constrain a model of the early 
Universe must satisfy. Present data constrain the non-gaussianity of the perturbations to be below 
|17j and this also needs to be explained. In the next section we extend our procedure to 
the non-linear case, so that perturbations can be followed to our present expanding phase also at 
non-linear level. If our assumption of neglecting the exponentially small terms is disregarded, the 
predictions for the higher moments of the perturbations will again depend on the unknown physics 
at the bounce. 

5 Non-linearly through the bounce 

Given our assumptions, we have shown in the last section that a bouncing model is not compatible 
with data. It seems therefore useless to proceed with our analysis at the non-linear level. We do it 
for several reasons. First of all, it clarifies the physical meaning of our argument, making it clear 
that it is not based on any linear approximation. Second it shows the analogy with the similar 
problem of following perturbations at the non-linear level in inflation, for example to calculate the 
higher moments of the spectrum. Finally it is interesting to understand how predictions can be 
made even if the cosmological evolution goes through a high curvature and strong coupling regime 
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at the bounce, which at first sight seems to jeopardize any perturbative approach. 

We have already stressed many times that after a mode leaves the horizon its effect becomes 
rapidly negligible for a local observer, analogously to what happens in inflation. A perturbation 
induces an intrinsic curvature of comoving surfaces. This effect blue-shifts quite slowly, as a~ 2 , 
so that it becomes rapidly irrelevant with respect to the unperturbed contraction: k 2 / (a 2 H 2 ) ~ 
t 2 ~ 2p . Once the curvature has become negligible, a local observer will see a homogeneous, flat but 
anisotropic Universe: perturbations leaving the horizon have left an anisotropic initial condition. In 
synchronous gauge, the metric of an anisotropic homogeneous, flat space can be put in the form 

ds 2 = -dt 2 + a 2 (t) e 2 ^ t) dx i dx i , ^ ft = , (37) 

i i 

where @i(t) describe the anisotropic expansion. Their evolution is simply given by 

A = Qa" 3 , (38) 
and their effect in the Friedmann equation blueshifts as a~ 6 , 

3M 2 H 2 = 3M 2 (~\ =p+ ^(c 2 + c 2 + c 2 ) . (39) 

These equations together with the equation of motion for the scalar field, <j) + 3H<j) + V'(4>) = 0, are 
all we need to describe the anisotropic evolution of the Universe. It is straightforward to check that 
eqs. (|32j) and (|33|) are a linearized solution of the equations above. The given non-linear description 
of an anisotropic Universe has been used in refs. j!21 118j to show that a contracting phase with 
w > 1 has an attractor towards which nearby solutions flow. All curvatures and anisotropics get 
diluted away, the fluctuations in the scalar drop to zero and we locally get back to the unperturbed 
solution up to exponentially small terms. In fact, even though it is not easy to get explicit solutions 
because the perturbation induces also fluctuations in the equation of state w, we see that the 
qualitative features we got in linear analysis remain valid. For instance in eq. (|39|) the unperturbed 
energy density goes as p ~ aT 2 ^ ', so that the anisotropic terms give a relative correction going like 
t 2 6p , which becomes irrelevant for p < 1/3. We therefore expect that the unperturbed 
solution is recovered up to corrections that drop off like powers of t 1_3p , in close analogy with the 
linear case. 

Following the same argument we used in the previous section we can neglect the exponentially 
small terms. Although distant regions of the Universe experience the same unperturbed history, as 
in the linear case their local metric is characterized by different, locally unobservable, integration 
constants. In general these constants describe the rescalings of the three spatial coordinates (along 
axes that can rotate from point to point) and thus the metric takes the form 

ds 2 = -dt 2 + e 2 f(*) e 2 7^) a 2 {t)dx l dx j , 7ii = . (40) 

The anisotropic piece 7^ comes from the fact that we have not restricted the analysis to scalar 
perturbations, as we instead did in the previous sections. Given a generic scalar-tensor initial 
configuration, going towards the bounce the system is driven to the form eq. (|40[) up to exponentially 
small terms. The prescription to follow the long wavelength perturbations to the expanding phase is 
now clear: as in the linear case, the metric will have the same form eq. (|40|) after the bounce with the 
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same space-dependent constants. They cannot change between the contracting and expanding phase 
because their variation would imply observable deviations from the unperturbed solution, while we 
know that every comoving observer is following the unperturbed history up to exponentially small 
terms. 

In the long- wavelength approximation, i.e. for modes well outside the horizon, the variables 
£ and 7 are non-linearly conserved independently of the phases the Universe goes through. This 
property is crucial to follow the perturbations created during inflation and to make the non-linear 
observables, for instance the 3-point function of density perturbations, independent of the unknown 
physics of reheating The same property holds also in our case, always assuming that the 

unknown bouncing phase is not sensitive to the exponentially small deviations from the unperturbed 
solution. Notice that the same assumption is implicitly made in the case of reheating. 

6 Conclusions 

The standard cosmological model has been extremely successful at explaining the evolution of the 
structures we observe in the Universe, from the early epochs we measure using the anisotropics in 
the cosmic microwave background to the large scale structures in the distribution of galaxies we 
see locally. However this model requires a mechanism to create initial perturbations that correlate 
points with separations larger than the horizon. Such initial seeds will then grow under the action 
of gravity to form the structures we observe. 

The natural assumption is to postulate that the initial seeds were created during an early phase 
in the history of the Universe and somehow "stretched" outside the horizon. The "stretching" 
of perturbations can be accomplished in two different ways. One can have a period of accelerated 
expansion in which the length scale of perturbations grows rapidly, becoming larger than the horizon 
which evolves only slowly. This is done in inflationary models. The other possibility is to have the 
Universe contract slowly so that the horizon shrinks faster than does the scale of perturbations. 

One fact is certain, whatever the scenario, accelerated expansion or slow contraction, this period 
has to end and give rise to the standard hot Big Bang phase. In inflationary scenarios the transition 
epoch is called reheating and could involve rather complicated physics. The beauty of inflation 
lays on the fact that the predictions of the model are independent of the details of reheating. The 
insensitivity to reheating stems from the fact that the wavelengths of interest to astronomy were 
extremely large compared to the horizon at the time and to the presence of a dynamical attractor. 
If the details of reheating were crucial to determining the predictions of inflation the scenario would 
not be that appealing. 

If perturbations were created during a contracting phase, with the Universe evolving towards a 
big crunch, the connection between the contracting phase and subsequent hot Big Bang phase may 
seem more problematic. The Universe has to go through a bounce with curvatures diverging as 
it approaches. This has been considered problematic by many detractors of such scenarios; in our 
mind however such criticisms are a bit unfair. Just as for inflationary models the important point 
is that their predictions are insensitive to the details of reheating, the real question for bouncing 
models is not how the bounce happened but whether or not one can show that their predictions are 
independent of the yet unknown UV physics responsible for the bounce. If one could argue that 
predictions are independent of the details of the bounce itself then bouncing models would be in 
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almost as good a shape as inflation. 

In this paper we have studied whether the predictions for density perturbations in the ekpy- 
rotic/cyclic scenario are robust, i.e. independent of the details of the unknown bouncing phase. The 
presence of a dynamical attractor in the contracting phase leads to a simple and natural assumption, 
namely that the physics of the bounce is independent of the exponentially small deviations from the 
unperturbed solution present as the bounce approaches. Even though seldom stated, this is the exact 
same assumption one makes in inflationary models about reheating. This assumption allowed us to 
calculate all the statistical properties of the fluctuations in the expanding phase. Unfortunately the 
spectrum is not scale invariant and thus incompatible with the data. 

If this assumption is evaded, all predictions (scale invariance of the spectrum, level of non- 
gaussianities, etc.) strongly depend on the physics in the high curvature regime of the bounce. In 
this case, predictions cannot be made without a full resolution of the singularity. In a sense our 
conclusions are more pessimistic than just the statement that in order to make bouncing models 
viable one needs to understand the details of the bounce. We have shown that from the perspective 
of a local observer the bounce has to be exponentially sensitive to the initial state. 

In the literature many prescriptions have been proposed to match the two independent modes in 
the contracting phase to those in the expanding one El E3 EH 120] • These prescriptions usually 
involve variables (like the Newtonian potential <!>) which diverge going towards the bounce and thus 
hide both the perturbativity of fluctuations and the existence of a dynamical attractor, features 
which are instead manifest in the synchronous gauge. Moreover, as different variables like $ and 
£ behave very differently in the contracting phase, one is lead to different "natural" prescriptions 
depending on which variable and gauge is used. We think that these mathematical prescriptions 
are not based on well motivated physical assumptions, and we consider the attractor as the only 
physical guide across the bounce. In fact it should be now clear that prescriptions giving a scale 
invariant spectrum describe a bouncing phase which is sensitive to exponentially small departures 
from the unperturbed solution, or in which all relevant modes come back into the horizon. In both 
cases results are anyway completely dependent on the unknown UV physics. 

For example some of the proponents of the ekpyrotic/cyclic scenario have developed a matching 
procedure based on analytic continuation in the full 5-d setup [U (see also ref. |22j). The result 
is that £ has a jump at the bounce proportional to the comoving energy density perturbation e m , 
AC ~ e m /(k p hL) 2 , where L is a UV length scale characterizing the physics of the bounce. The 
variable e m goes to zero during the contracting phase as t 1 ~ 3p and it is therefore exponentially small 
at the bounce. Therefore from the 4d perspective this prescription describes a bounce which is 
exponentially sensitive to the incoming state; e m is in fact multiplied by a huge ^-dependent term. 
This extreme sensitivity to initial conditions is clearly problematic; for instance it amplifies also all 
preexisting anisotropics which were diluted away during the contracting phase thanks to the presence 
of the attractor. 

We conclude stressing that the situation is slightly different for pre-Big Bang scenarios [3]. 
In this case there is no potential for the scalar (dilaton) and only its kinetic energy is driving 
the (Einstein frame) contracting phase; this corresponds to the limit p — > 1/3 in our solutions 
of section |2 The kinetic energy blue-shifts as a -6 , exactly like the contribution of anisotropies 
to the Friedmann equation (|.39|) . This implies that we are somewhat borderline with respect to 
the existence of the attractor. In fact the anisotropic perturbation in synchronous gauge and the 
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fluctuation in the scalar field, instead of going to zero, diverge logarithmically. The situation now 
resembles an inflationary cosmology in the presence of isocurvature perturbations; separated regions 
of the Universe are characterized by a different mixture of two independent modes and therefore 
go through the bounce in a different way depending on the relative size of the two solutions. Our 
conclusions cannot be applied and it is necessary to specify the details of the high curvature regime 
to evolve the fluctuations across the bounce. 
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Appendix 

A The quadratic action for scalar fluctuations 

The complete action for gravity and a minimally coupled scalar field <j) is 

S = i J d 4 xy^ [R - {d^f - 2V(cj>)] . (41) 

(For notational simplicity we are setting Mp = 1; at the end of the computation Mp can be restored 
by dimensional analysis) . In order to find an action for the scalar fluctuations of this coupled system 
around a given background it is convenient to work in the ADM formalism. This has been done 
in ref. ^3], and in this section we briefly review the results derived there. In terms of the ADM 
variables the metric is 

ds 2 = -N 2 dt 2 + h ij (dx i + N i dt)(dx j + N j dt) , (42) 
and the action eq. (|41|) reads 

NRW - 2NV(<t>) + N~ 1 (E ij E lj - E 2 ) + N~ x (</> - N l d l( f) 2 - Nh ij drfdjcj) . 

(43) 

Here R^ is the intrinsic curvature of spatial slices, while Eij is related to the extrinsic curvature, 

Eij = \ {hij - V t Nj - VjNi) , Kij = N- l Eij , (44) 

and E = E l {. From the action above it is clear that N and N l are just Lagrange multipliers. One 
can then solve for them the constraint equations, express them in terms of the other degrees of 
freedom, and substitute their value back in the action. At this stage it is useful to pick a gauge. We 
choose the (-gauge, defined at linear order in the fluctuations by 

</>@,t) = Mt) , hij = a 2 (t)(l + 2C)<% , (45) 



S = I / d 4 xVh 
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where 4>o(t) an d a(t) are the background solutions for the scalar and the scale factor, and we have 
neglected tensor modes. In this gauge the constraint equations read 

ViiN-^E'j-SiE)) = 0, (46) 
-2V{4> Q )-N- 2 {E ij E i ^ -E 2 ) -4>l = 0, (47) 

whose solution at linear order is 

Plugging this into the action above and expanding up to second order one gets a quadratic action for 
the only scalar degree of freedom (. After integrating by parts and using the background equations 
of motion {i.e. the Friedmann equation and the equation for the scalar) one finally gets 

S=\J dtd 3 x || [a 3 ( 2 - aidiCf] =-\J d A x^ j| WO , (49) 

where is the background FRW metric. For the scaling solution discussed in the text the quantity 
4> 2 /H 2 is a constant, and the action above reduces simply to eq. (|T2*|) . 

In the alternative gauge eq. (|7j). which keeps 5<p as the scalar dynamical variable, the result for 
the action is not as simple as in the (-gauge. This is because the relation between ( and 5(j> is in 
general time-dependent, ( = — (H/4> ) 8(j), so that time derivatives acting on ( in the action above 
also produce time derivatives of the factor H/<pQ. However, for the scaling solution we are interested 
in such a factor is time-independent, so that again the quadratic action takes the very simple form 
eq. (JBJ). 



B From (-gauge to synchronous gauge 

In this section we explicitly derive the gauge transformation that goes from ("-gauge to synchronous 
gauge. In (-gauge the metric for a general scalar fluctuation is at linear order (see the previous 
section) 



9oo 



H 



90i 
9ij 



-0, 



H 



it 



-di. 



(50) 

(51) 
(52) 



2H 2 8 2 
(1 + 2C)<% • 

Under a gauge transformation of parameter £ M the metric transforms as g^ v — > g^ u + + V;,£ M , 
where the V's are the covariant derivatives associated to the unperturbed FRW metric. The non-zero 
Christoffel symbols in a FRW geometry are 



A) 



a 2 H5i 



H5 1 : 



so that the explicit transformation of the metric components is 



500 — > 9oo + 2£o , 

9oi -> 9oi + 4 + - = g 0i + a 2 — -| + d^ , 

dt a A 

9ij -> 9ij + diij + djii - 2a 2 H5ijio 



(53) 

(54) 
(55) 
(56) 
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We want to end up in synchronous gauge, which is defined by goo = — 1 and go % = 0. The explicit 
form of goo combined with its transformation law tells us that we must choose 



to 



H 



(57) 



To be completely general we should also add to to a generic function of x, but it will turn out that 
it is not necessary, so we set it to zero. 

Then we want to set to zero the space-time components goi- In order to do so must satisfy 



dt a 2 a 2 



so that 



ii = a 2 di / dt 



j 1 



j-dt" 

o H 



(58) 



(59) 



The last step is to combine eqs. (|521 156(1 and read off the spatial metric in synchronous gauge, 



gij -+ a 2 { ( 1 + 2C - 2H j jjdt'^J 5ij + 2didj / dt' 



1 



2H 2 8 2 ^ a 2 \H 



1 c 



, (60) 



where we used the definition of H, eq. (|51|1. Specializing to our scaling background a oc t p , H = p/t, 
(fro = y/2p/t we get precisely eq. (|32|) . 

Finally we want to express the scalar <fi(x, t) in synchronous gauge. Under a gauge transformation 

<f>{x) -> 0(jc) + d^ix) ■ e . (61) 
Since in ^-gauge the scalar is unperturbed, 4>(x,t) = (j>o(t), in synchronous gauge we simply have 

ft £ 



= <j)(x, t) - (j)o(t) = 00 e° = -00 CO = -00 f ±<U' 

JO 



H 



(62) 



For our scaling solution this reduces to eq. 
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